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ABSTRACT 


Calculations  are  made  to  estimate,  as  a  function  of  distance,  the  radiation 
dose  levels  which  would  occur  during  full  capacity  operation  of  the  Kray  Pulse 
Reactor  Assembly  (APRA)  scheduled  for  construction  at  Aberdeen  Proving  Ground, 
Maryland.  Calculations  of  the  neutron  dose  are  made  by  two  methods.  The  first 
calculation  (Method  A)  is  made  by  utilising  an  experimentally  determined 
effective  dose  wan  free  path,  which  probably  Includes  the  effect  of  ground 
scattering.  The  second  calculation  (Method  B)  is  made  by  utilizing  the  result 
of  a  Monte  Carlo  calculation  of  neutron  transmission  in  air.  The  results  of 
this  calculation  indicate  considerably  higher  dose  levels  than  indicated  by 
Method  A.  Experimental  data  are  available  which  substantiate  the  result  of 
Method  A.  The  gamma  dose  at  distances  of  interest  Is  found  to  be  unimportant 
when  considered  in  light  of  the  uncertainty  in  the  neutron  dose. 
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I.  INTRODUCTION 


In  order  to  determine  the  degree  to  which  operation  of  the  Arny  Pulse 
Reactor  Assembly  (AFRA)  at  the  proposed  Aberdeen  Proving  Ground  site  will 
disrupt  other  nearby  activities,  it  is  necessary  to  calculate  the  radiation 
dose  which  will  be  delivered  at  various  distances  from  the  nuclear  reactor. 

In  the  case  of  an  unshielded  reactor  such  as  AFRA,  the  exclusion  area 
will  be  determined  by  the  dose  produced  by  radiation  emanating  from  the  reactor. 
The  proposed  site  is  in  the  industrial  area  of  the  Proving  Ground.  Access  to 
this  area  by  the  public  and  by  all  Aberdeen  Proving  Ground  employees  except 
those  who  have  industrial  area  passes  is  restricted.  The  following  proposals  have 
been  to  further  restrict  access  to  the  reactor  site  and,  at  the  same  time, 
to  provide  adequate  protection  to  insure  the  safety  of  persons  not  directly 
under  the  administrative  control  of  the  reactor  staff. 

a.  There  should  be  a  barbed  wire  fence  with  appropriate  signs  attached 

at  a  distance  such  that  the  maximum  weekly  dose  outside  of  this  fence  is  10  mrem. 

b.  TCiere  should  be  a  security  fence  (6  or  8  foot  fence  with  barbed  wire 
barrier  on  top)  at  a  distance  such  that  the  maximum  daily  dose  outside  of  this 
fence  is  10  rem. 

Members  of  the  Hazards  Evaluation  Branch  of  the  AEC  have  indicated  that 
these  restrictions  would  probably  be  considered  adequate  for  exclusion  of  "non- 
occupational"  personnel. 

II.  CALCULATION  OF  NEUTRON  DOSE 

'Hie  Godlva  II  neutron*  leakage  spectrum  is  used  as  a  source  model  for  the 
calculation  by  Method  B.  See  Reference  1.  This  spectrum  is  Included  as 
Figure  1.  The  total  source  is  assumed  to  be  80  x  10^7  fissions,  and  the  number 
of  leakage  neutrons  per  fiBslon  is  assumed  to  be  1.4  (Reference  2,  p.  9) •  This 
source  is  equivalent  to  40  pulses  per  week  at  2  x  10^7  fissions  per  pulse. 

2  x  1017  fissions  is  the  probable  maximum  yield  which  the  AFRA  will  deliver 
in  one  pulse.  The  dose -distance  calculation  has  been  made  by  two  different 
methods. 

*4he  7  dose  is  estimated  at  about  2  per  cent  of  the  neutron  dose  at  distances  of 
interest.  See  Section  III. 
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EXTRAPOLATED 
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NEUTRON  ENCROY  (MEV) 


A.  Method  A 


A  simplified  recipe  for  calculating  the  dose  in  air  at  several  mean 
free  paths  from  a  neutron  source  has  been  proposed  (Reference  3).  The  proposed 
method  is  as  follows:  (l)  Use  an  infinite  mean  free  path  (no  absorption)  for 
the  first  300  feet.  (2)  Use  a  mean  free  path  of  212  yards  from  300  feet  on  out. 
Changes  in  air  density  would  cause,  presumably,  a  change  in  mean  free  path. 

A  calculation  was  made  using  the  above  method  with  an  infinite  mean  free 
path  for  the  first  100  meters.  This  method  predicts  10  mrem/wk  at  1110 
meters  (1215  yards). 

B.  Method  B 

A  more  complex  calculation  was  made  using  the  data  presented  by  Mehl 
(Reference  4).  Mehl  gives  dq^e -distance  data  in  air  for  nine  monoenergetic 
sources.  The  Godiva  spectrum  was  approximated  by  building  up  a  spectrum  composed  of 
the  nine  monoenergetic  soirees  given  by  Mehl.  An  appropriate  number  of  neutrons 
was  assigned  to  each  source  according  to  the  number  of  neutrons  of  the  Godiva 
spectrum  in  the  energy  range  corresponding  to  monoenergetic  source.  For  each 
of  the  nine  monoenergetic  sources,  Reference  4  gives  the  number  of  neutrons 
of  a  given  energy  at  various  distances  from  the  source.  For  Instance  if  one  is 
interested  in  finding  the  dose  received  at  1000  meters  from  a  source  of  1  Mev 
neutrons,  data  are  given  to  enable  one  to  determine  the  number  of  neutrons  in 
each  of  several  energy  ranges,  ltoe  neutron  dose  at  a  distance  R  is  determined 
by  accumulating  the  number  of  neutrons  arriving  at  R  in  each  energy  range  from 
all  source  energies.  The  number  of  neutrons  of  a  given  energy  required  to 
make  one  rem  is  given  in  Reference  5*  These  data  are  presented  in  Figure  2. 

A  sample  calculation  (with  R  =  1500  meters)  is  given  on  pages  6  and  7.  The 
calculation  is  based  on  a  toted  source  of  one  neutron.  Source  neutron 
energies  are  distributed  as  indicated  in  Figure  1.  Symbols  used  are  defined 
as  follows: 

i  -  subscript  used  to  denote  source  neutron  energy. 

AEj  -  Energy  range  of  the  source  neutrons 

F.  -  Fraction  of  the  source  neutrons  in  energy  range  AE.  as  determined 
from  Figure  1. 
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NEUTRON  ENERGY  (MEV) 


FIG.  2.  n/cm2  >  IO*7 


R 

J 

AE 


J 

ij 


'J 


-  Range  (distance)  from  the  source  in  cm. 

-  Subscript  used  to  denote  neutron  energy  range  at  R. 

-  Energy  range  of  the  neutrons  at  R 

-  Number  of  neutrons  per  square  cm  at  R  in  energy  range  AEj  which 
started  with  a  source  energy  in  AE^^  (assuming  total  source  of  one 
neutron) 

p 

-  Conversion  factor  for  determining  the  number  of  neutrons/cm  of 
energies  AEj  necessary  to  make  one  rem. 

>  Neutron  dose  at  R  due  to  neutrons  with  energies  in  energy  range  AE.. 

U 


The  air  density  used  by  Mehl  is  1.07  grams/liter.  The  mean  July  air  conditions 
for  AFQ  are  as  follows: 


T 

avg 

B 


24°  C 

761.5  mm  Hg 


Rel.  Hum  =  71  per  cent 

These  conditions  indicate  a  mean  July  (when  the  monthly  mean  density  is 
lowest)  density  of  1,101  grams/liter.  A  density  correction  1b  applied  according 
to  the  following  formulas: 


V  PlV  *  (  V  2)2  D2(  2R2) 

where  p  =»  air  density 

R  a  distance  from  the  source 
D  *  dose 

•PI 

In  the  calculation  above,  we  have  R^  =  1500  meters,  ■  2.42  x  10  rem/source 

neutron,  P  1  -  I.07  grams/liter,  P  0  »  1.101  grams/liter.  Using  these  values  in 
the  above  equations,  one  gets:  Rg  -  1559  meters,  Dg  ■  2.95  *  10  rtnq/ source 

neutron  (1.4  source  neutrons/f lesion)  a  4. 15  x  10-21  rem/ fission.  Additional 
dose-distance  data  are  shown  in  tabular  form  in  Table  II  and  in  graphical 
form  in  Figure  5.  The  "Inverse  square  loss  only"  curve  is  shown  for  comparison 
purposes  only. 
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NEUTRON  DOSE  (REM/FISSION  ) 


TABLE  I  (Continued) 


z*« 

1  3 

/_  / _ c 


*  2 
(n/ca  i 

WKm^fR, 

xlO15 

xlO-7 

xlO22 

42.17 

79.0 

.534 

79-75 

53.0 

2.417 

15.65 

7.8 

2.032 

29.U 

4.0 

7.276 

12.4? 

2,74 

4.544 

6.70 

2.64 

2.538 

10.63 

2.84 

3.613 

2.366 

2.76 

.858 

.506 

2.6 

-^2L 

Total  doee/n 


24.2  z  lO"82  rw/n 


C.  Tabulation  of  results 


TABLE  II 

Dose-Distance  Data  for  Neutrons 


Calc.  Method 

Range 

(Meiers) 

Dose 

(rem/f  lesion) 

A 

960 

-21 

3.75  x  10  " 

A 

550 

9*41  x  10’20 

A 

200 

4.31  x  10‘l8 

B 

1359 

4.13  x  10“ 21 

B 

906 

6.83  x  10"20 

B 

272 

7.00  x  10"18 

C* 

1200 

2.02  x  10”19 

C 

900 

3.58  x  10" 19 

c 

500 

1.16  x  10“l8 

c 

300 

3.23  x  10_l8 

c 

180 

8.96  x  10"18 

C(and  A) 

100 

2.90  x  10"17 

In  addition  to  the  above  data  there  are  bosks  experlaental  data  available 
from  the  Sandla  Pulse  Reactor.  The  following  data  are  reported  (Reference  14) 
for  the  reactor  operated  at  the  outdoor  site. 


Dose 

(rea/fission) 


250  8  x  10"17 

1750  5  x  10" 20 

These  data  agree  very  closely  with  calculation  A  as  shown  in  Figure^. 


♦Calculated  dose  assuming  no  attenuation  other  than  inverse  square  loss. 
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III.  CALCULATION  OF  GAtftiA  DOSE 


A.  Fission  Gammas 

The  total  fission  gamma  source  (delayed  plus  prompt)  Is  assumed  to  be 
13.1  Mev/fission.  See  Reference  6,  p.  336.  The  distribution  of  gammas  is 

taken  as  n(E)  =  A  exp  where  c  =  1.16  after  Reference  7.  By  setting 

2 

A  =  c  ,  one  normalizes  the  distribution  to  a  total  source  of  1  Mev.  Note 
that  the  average  energy  is 


If  c  =  1.16,  E  =  0.86  Mev.  For  prompt  gammas,  c  is  sometimes  taken  as  —  1.4 
(e.g. ,  Reference  7;  also  the  estimated  slope  of  the  straight  line  drawn  through 
data  of  Figure  2,  Reference  6  is  such  that  c  —  1.4).  This  gives  E  for  prompt 
gammas  of  ~  .7  Mev.  However,  Reference  6,  p.  369  gives,  for  prompt  in 

the  energy  range  .3  to  10  Mev,  7*2  +  0.8  Mev/fission  carried  by  7*4  +  0.8 
photons/fission,  or  an  1  of  2  1.0.  By  manipulation  of  data  presented  in  Reference 
6,  one  finds  that  E^—  0.95  Mev  averaged  over  fission  gammas  emitted  promptly, 
short  delayed  and  long  time  delayed.  If  E^  —  0.95  Mev,  the  calculations  given  here 
indicate  a  gamma  dose  which  is  too  small  by  5  to  10  percent  when  one  takes  into 
account  the  fact  that  the  nitrogen  capture  gamma  dose  is  dominant  at  large  distances. 
Since  it  will  be  shown  that  the  gamma  dose  is  small  in  comparison  with  the  neutron 
dose  for  the  situation  under  consideration,  no  further  attempt  will  be  made  to 
resolve  this  apparent  contradiction. 

2  _cE 

Using  the  distribution  n(E)  =  13.1  c  exp.  ,  the  energy  range  0  to  8  Mev 
is  broken  down  into  energy  groups  of  l/2  Mev  width  as  shown  in  Table  III.  Bxe 
number  of  gammas  in  each  group,  niQ,  is  assigned  the  median  energy  of  the  energy 
range,  E.. 

*  • 


n^Q  -  n(E)dE  =  13.1  c 

The  number  of  photons  in  each  energy  group 
to  the  formula 


£exp.  fCEl)  -  exp  (“CE2)J 
(nio)  is  then  attenuated  according 


-  QioBgC^ZjE)  exp(’^sZ) 
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where  Bj,(pflZ,E)  is  the  energy  buildup  factor  for  a  plane  monodirectional  gamma 

source  in  uranium  obtained  from  Reference  8,  p.  146.  n^f.P. )  is  the  number 

of  photons  due  to  fission  products  in  energy  group  i  which  escape  from  the 

reactor,  u  (E)Z  is  the  "effective"  number  of  relaxation  lengths  for  the 
8 

photon  energy  E.  This  factor  takes  into  account  the  fact  that  many  of  the 

photons  produced  are  absorbed  in  the  reactor  core  itself.  The  relationship 

between  the  number  of  source  radius  relaxation  lengths,  p  (E)d,  and  the 

s 

"effective"  number  of  relaxation  lengths  is  given  in  Reference  9,  p.  75* 

Hie  radius  of  the  reactor  core,  d,  is  assumed  to  be  4.5  inches  of  which  7 
percent  (of  the  radius)  is  gaseous  coolant  channel  which  we  will  assume  to  be 
non-attenuating.  These  values  were  estimated  from  data  given  in  Reference  10, 
p.  42.  Hie  reactor  core  is  90  percent  U  and  10  percent  by  weight  Mo.  Hie 
linear  absorption  coefficients  are  calculated  by 


ufi(E)  =  0.95  To.9  +o.i  ^  PuMq 
where  L  J 

*>  the  mass  absorption  coefficient  for  U, 

^n*to  a  the  1,1688  atsorPtiOQ  coefficient  for  Mo, 

Pt#to'  the  denslty  t*ie  U-Mo,  is  taken  to  be  17.08  gm/cm^  after 

Reference  11,  p.  4.  Values  of  the  mass  absorption  coefficients  u  (E)  are  taken 

SI 

from  Reference  12.  Thus  the  number  of  photons  in  each  energy  group  escaping 
fro--  the  reactor  1b  determined. 


Nitrogen  Capture  Gammas 

B.  This  source  is  directly  dependent  on  the  number  of  leakage  neutrons  per 
fission,  T|.  The  value  of  1}  is  taken  as  1.4  neutrons/fission  (Reference  2,  p.  9)» 
Hie  number  of  Mev  of  nitrogen  capture  (N.C. )  gammas  produced  per  fission,  U^, 
is  computed  as  follows: 


U, 


N 


NC 


where 


sair  / 


air 


rT^c 


«  ratio  of  neutron  absorption  in  nitrogen  to  that  in  air. 


air  /  abs 


>  ratio  of  radiative  cature  (N.C.)  to  absorption  in  nitrogen. 
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y, j,  ■  average  number  of  photons/N.C. 

Eq  »  average  energy  of  photons  emitted  in  H.C. 
The  values  used  In  this  calculation  are: 


1  ) 

°air  / 

-) 

i/N 


Zz 

0aj 


Eo  = 


=  1.0 


.1 

lTBT- 


m  -  1.44 
=  6.44 


5^2 

48.46 


U, 


NC 


7-53 

0.71  Hev/fission 


The  distribution  of  gammas  emitted  In  neutron  captures  In  nitrogen  Is  given 
in  Table  IV  (taken  from  Reference  j). 

TABLE  IV 

Photon  Emission  Energies  and  Probabilities  of  Heutron  Capture  by  Iltrogen 


k 

Energy 

\>  MeV 

Relative  Ho. 

*k 

P 

k 

1 

1Q.816 

1.0 

.01456 

2 

9.156 

.09 

.00151 

3 

8.278 

.19 

.00277 

4 

7.356 

.56 

.00815 

5 

7.164 

•19 

.00277 

6 

6.518 

•90 

.01510 

7 

5*554 

1.5 

.02184 

8 

5.287 

2.3 

.05349 

9 

4.485 

.8 

.01165 

Total 

7-53 

9 


pk  * 


^(•73.) 


No.  of  photons  of  Ener^1  per  fission 


The  number  of  photons,  P^,  of  energr  E^  is  then  assigned  to  the  median  energy 
of  the  particular  energy  group  in  which  the  energy  E^  fails  as  indicated 
in  Table  III. 


The  nitrogen  capture  gammas  are  treated  as  a  point  source  positioned  at 
the  reactor  center.  This  approximation  is  Justified  on  the  basis  that  in 
this  calculation  we  are  only  interested  in  the  dose  at  1000  yards  or  more  and 
the  great  majority  of  all  neutrons  are  absorbed  well  within  this  distance. 


C.  Composite  Gamma  Dose 

The  number  of  gammas  in  each  of  19  ener^  ranges  between  0-11  Mev, 

n^( corap)  =  n  f  is  determined  by  adding  the  contributions  from  fisBion  gammas 

and  N.C.  gammas.  The  values  of  n^  are  listed  in  Table  III.  Also  listed  in 

Table  III  are  values  of  r(x,E. )  =  r,  which  is  the  energy  in  Mev  deposited  per 
x  11 

cm  in  air  at  a  distance  x  in  yards  from  a  source  of  one  photon  of  energy  E^. 
Values  of  r^  are  given  in  Reference  7  and  are  calculated  by  the  formula 


r 


i 


exp(-u^x) 

4*x2 


where  Brl  is  the  dose  buildup  factor  in  air  for  photons  of  energy  E^,  p*  is  the 
linear  absorption  coefficient  not  including  corapton  scattering  component,  ^ 
is  the  total  linear  absorption  coefficient  including  the  compton  scattering 
component,  and  x  is  the  distance  from  the  source  in  cm.  The  buildup  factors 
are  taken  from  Reference  7  as  before.  The  air  density  for  the  gamma 
calculations  is  taken  as  1.205  gmA  rather  than  l.lSl  gms/l  since  data  are 
already  available  at  1.205  gm/l.  This  small  difference  will  result  in  very 
little  inaccuracy  of  the  final  results. 


and 


Prom  Table  III,  we  see  that 

V  n^ilOOO,  E1)  =  3.13  x  10 

n^  (1500, Ei)  =  3.79  x  10 


^Mev/cm^  -  fission 
Mev/ cm^-f ission 
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Using  the  conversion  factors, 

(l.6  x  10‘6  ergs/Mev)  (  (hIBIx  IO^  “1*545  X  lO’^/Mev/c*5, 

one  obtains  the  gamma  dose  per  fission  shown  in  Figure  4.  For  a  2  x  lO1^ 
fission  pulse  the  gamma  dose  is  ~0.1  mr  at  1000  yds  and  0.01  mr  at  1500  yards. 

IV.  ANALYSIS  OF  CALCULATION 

A.  Results 

The  gamma  dose  is  less  than  2  per  cent  of  the  neutron  dose  (as  calculated 
by  Method  B)  for  distances  out  to  about  1500  meters.  Therefore,  we  may  neglect 
the  gamma  dose. 

The  distance  at  which  10  mrem/wk  will  be  accumulated  from  a  weekly  source 
IT 

of  80  x  10  '  fissions  is  1215  yards  based  on  calculation  Method  A.  The  corresponding 
figure  for  Method  B  is  1700  yards. 

The  distance  at  which  10  ren^ay  will  be  accumulated  is  195  yards  by  Method 
A  and  315  yards  by  Method  B. 

B.  Factors  Bearing  on  the  Results 

It  is  believed  that  the  results  of  Method  A  are  more  correct  than  the 
results  of  Method  B.  This  contention  is  supported  mainly  by  the  experimental 
data  presented.  Factors  having  an  effect  on  either  or  both  of  the  two  methods 
will  be  enumerated.  The  applicable  calculations!  method  Is  indicated  by  the 
letter  (s)  preceding. 

A,B  1.  The  radiation  leakage  from  the  AFRA  may  be  more  (or  less)  per  fission 
than  that  for  Godlva  (which  was  used  in  these  calculations). 

B  2.  Thermal  and  resonance  neutrons  arriving  at  the  detector  below  200  ev 
energy  have  not  been  accounted  for  in  this  calculation. 

A,B  3*  Ihe  AFRA  neutron  leakage  spectrum  may  not  be  the  same  as  the  Godlva 
spectrum  (which  was  used  in  these  calculations). 

B  4.  This  calculation  is  made  for  an  infinite  medium  of  air  with  a  point 

source  of  neutrons  (and  gammas).  The  ground-air  Interface  may  decrease 
the  dose  at  large  distances  from  the  reactor  particularly  when  the 
reactor  is  operated  close  to  the  ground.  Data  are  available  for  games 
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ROENTGENS/FISSION 


(Reference  13,  p.  89);  however,  no  data  has  been  found  for  neutron 
transmission.  The  non-nitrogen  capture  gamma  dose  will  probably  be 
reduced  by  about  a  factor  of  2  or  3  at  a  distance  of  1500  meters.  It 
would  not  necessarily  follow  that  this  same  figure  can  be  used  for 
neutrons,  however. 

A,B  5*  It  may  not  be  necessary  to  pulse  the  reactor  once  per  hour. 
Nevertheless,  it  would  be  desirable  to  be  able  to  do  so. 

A,B  6.  The  APRA  may  not  be  capable  of  producing  a  pulse  which  is  13  times 

as  large  as  the  nominal  Godiva  II  pulse  of  1.5  x  10^  fissions/pulse. 

B  7.  Excessive  moisture  in  the  air  will  reduce  the  dose  level  below  that 
calculated  here.  Mehl's  calculations  assume  a  constant  air  density 
of  1.07  x  10 gna/cB?  with  relative  atom  fractions  of  nitrogen,  oxygen, 
and  hydrogen  of  O.7778,  0.2135  and.  O.OO85  respectively.  APG  air 
contains  more  than  twice  the  amount  of  moisture  as  that  used  in  Mehl's 
calculations  (assuming  mean  July  conditions  for  APG,  the  atom  fraction 
for  hydrogen  is  .0209). 

B  8.  The  data  used  to  calculate  the  doBe  by  Method  B  are  estimated  to  be 

accurate  to  within  +  50  per  cent  according  to  the  author  of  Reference  h 
(assuming  an  infinite  medium  of  air). 

A,B  9*  Any  attenuation  of  neutrons  in  the  walls  of  the  reactor  building  or 
experimental  equipment  positioned  around  the  reactor  will  reduce  the 
dose  at  distances  of  interest. 

Of  the  above  items,  perhaps  only  the  last  one  may  be  expected  to  result 

In  a  significant  change  in  the  predicted  dose. 


FREDERICK  H.  GREGORY  ' 
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